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Abstract 



o 
o 

(N 

P5 . The vacuum solutions in brane gravity differ from those in AD by a number of additional terms 

^ ' and reduce to the familiar Schwarzschild metric at small distances. We study the possible roles 

that such terms may play in the precession of planetary orbits, bending of light, radar retardation 
\l ' and the anomaly in mean motion of test bodies. Using the available data from Solar System 

experiments, we determine the range of the free parameters associated with the linear term in the 
metric. The best results come from the anomalies in the mean motion of planets. Such studies 
O , should shed some light on the origin of dark energy via the solar system tests. 

in . 

bJO, 
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1 Introduction 

> 

\ In the past two decades, various attempts have been made to understand the nature of the so-called 

■ dark energy and dark matter which would be required to explain many observational data. For a 

review, see for example [1, 2, 3, 4]. Deviation of the galactic rotation curves from Newtonian gravity 
\ which occurs at distances larger than the Solar System scales [5], the velocity of galaxies in clusters 

' and the bending of light rays from galaxies and clusters [6] are the best evidence for the existence 

. of copious amounts of dark matter. Apart from the efforts made by a number of authors to account 

for dark matter by considering such elusive objects as massive neutrinos, axions and the Weakly 
Interacting Massive Particles (WIMPS) which is predicted by supersymmetric theories, the nature 
of dark matter is still unresolved. An interesting attempt to explain the galactic rotation curves was 
5^ \ made in [7] where the author considers dark matter as a galactic phenomena and proposes what is now 

known as the Modified Newtonian Dynamics (MOND) at low accelerations. However, since MOND 
is a non-relativistic theory, it cannot predict any relativistic phenomena such as the bending of light, 
the precession of planetary orbits etc. Solar System experiments have been increasingly relied upon 
to investigate the integrity of the foundation of general relativity in recent years, see for example 
[8]-[18]. 

In the recent past, models incorporating extra dimensions have become the focus of attention 
for investigating the nature of problems mentioned above. In these theories one considers a four 
dimensional world (brane) embedded in a higher dimensional manifold (bulk) through which only 
gravity can propagate. Ordinary matter is confined to the brane and cannot propagate through the 
bulk. The confinement is achieved, through the imposition of Z2 symmetry and use of the Israel 
junction conditions which relates the extrinsic curvature of the brane to the energy momentum of 
the matter. This method has predominantly been used in theories with one extra dimension. If the 
number of extra dimensions exceeds one, no reliable method for confining matter to the brane exists. 
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This is so since the requirement to define junction conditions is the existence of a boundary (brane) 
which cannot be defined if the number of extra dimensions is more than one. For example, a boundary 
surface in a 3D space is a surface with one less dimension whereas a line in the same space cannot be 
considered as its boundary. On the back of such concerns, model theories have been proposed where 
matter is confined to the brane through the action of a confining potential, without the use of any 
junction condition or Z2 symmetry [19]. In [20] the authors used the confining potential approach to 
study a brane-world embedded in a m-dimensional bulk. The field equations obtained on the brane 
contained an extra term which was identified with the X-cold dark matter. The same methodology 
was used in [21] to find the spherically symmetric vacuum solutions of the field equations on the brane. 
These solutions were shown to account for the accelerated expansion of the universe and offered an 
explanation for the galaxy rotation curves. 

In this paper, we focus attention on the consequences of the spherically symmetric vacuum so- 
lutions mentioned above when considering such questions as the precession of planetary orbits, the 
defiection of light rays in the Solar System, the time delay of signals in the Solar System and the 
mean motion of test bodies. In doing so, we obtain constraints on the free parameters appearing in 
the metric. This should hep us to accounting for the origin of dark energy via Solar System tests. 



2 The model 

Let us start by presenting the model used in our calculations. We only state the results and refer the 
reader to [20, 21] for a detailed derivation of these results. 

As was mentioned in the introduction, the brane-world model we invoke here differs from the 
usual Randall-Sundrum type in that no junction conditions or Z2 symmetry is used. One thus starts 
with the usual setup in which a AD brane is embedded in a 5 or, in general, n-dimensional bulk. 
Assuming that the brane is devoid of matter with no cosmological constant and that the bulk space 
has constant curvature, one arrives at the following equations [21, 22] 

G^i/ = Q fill J (1) 

where 

Q^, = [K%Kp, - KK^,) - ^ [K^pK^P - K^) 5^.. (2) 
We note that is an independently conserved quantity, that is 

Q'-^ = 0, (3) 

so that equation (1) satisfies the covariant conservation law. Equation (1) is the starting point from 
which a class of solutions were found in [22], representing a black hole. Thus, starting with the metric 

the following solution satisfies equations (1) and (3) 

e'' = e-'' = l---a\'^-2jr-l3^. (5) 
r 

This solution represents a black hole which we shall consider in the next section and use to analyze 
the behavior of a test particle in such a space-time. 
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3 The geodesic equations of motion 



The induced line element of the vacuum 4D space-time is 



where 



ds^ = A{r)dt' 



2m 



Air) 



r^dVL\ 



A(r) = l a>^-27r-/3^ 

r 



(6) 



(7) 



with 7 = a/3, m is the mass of the central object, a and (5 are constants and we have adopted the 
rclativistic units, c = G = 1. In this space-time, the geodesic equations of motion for a test particle 
are the Euler equations resulting from the lagrangian 
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A{r)t^ 



A{r) 



[e"^ + sin^ e^fP'^ 



(8) 



where a dot denotes differentiation with respect to the affinc parameter. Without loss of generality, 
let us consider the equatorial plane, ^ = ^, for a test particle. The lagrangian then becomes 



L = i 
2 



A(r)t 



12 



Air) 



2 -2 



Now, using the Euler-Lagrange equations, we obtain 

E 



A(r)' 



(9) 



(10) 



J 



ill) 



with J and E being constants. Along the orbit for a test particle 2L = 1 and for photons 2L = 0. 
Now, using equations (10), (11) in equation (9) we obtain 



1 



A{r) 



E^-r'' 



J2 



(12) 



where C = I and C = represent the massive and massless particles respectively. 

4 Precession of planetary orbits 

To find equation of motion of a massive particle in this space-time we use equation (12) (C = 1) and 
obtain 



.'2 



+ {l-P') u' = 



E^ + p'^-l 2mu „ „ s 27 2 



(13) 



where u = ^ and a prime denotes differentiation with respect to ip. Now by differentiating equation 
(13) with respect to (p we find the differential equation of the motion for a massive particle as 



a 



(14) 



dip"^ V / J2 J2y3 J2y2 

In order to solve this equation we consider Smu^, -^p^ and as perturbative terms, since these 
terms are much smaller than j^. This can be seen hy noting that, for example, "dmu^ j [m j J'^) = 



3 



(3/r^)(r^(/j)^ ?){rdipld£f l(? k, "iv^ jc? k, 7.7 x 10 ^ for Mercury. To first order then, this equation 
has a solution of the form 



u ~ 



1 + e cos 1 — 



1 



where 



(1-/^0 



1/2 



(15) 



(16) 



7+7. 

and e is the eccentricity. From equation (15) we find that r is a periodic function of ip with period 

1/2 - o„. „ XTl 



27r 



(l_/32)l/2 



3a2p4 ^p3 3^ 



> 27r. 



17) 



Since /?^,7 ^ 1, to first order approximation, the perihehon anomaly At/? of the particle after one 
revolution is found to be 



3m2 /J4 \ 



3j6a2 
2m4 



1 9m^ 



2J2 



(18) 



For 7 > 0, we find that the perihelion anomaly will increase compared to what one has for the 
Schwarzschild space-time [23, 24, 33]. In the case 7 < 0, depending on the parameters representing 
the particle, the perihelion anomaly may increase or decreases compared to what is predicted by the 
Schwarzschild metric. 



5 Deflection of light rays 

In this section we investigate the deflection of light in a space-time represented by the metric above 
in two cases, 7 < and 7 > 0. We assume that both the observer and the light source are far 
from the compact object whose presence would deflect the light rays. The light emitted from the 
source approaches the compact object, reaching the observer at infinity. Using equation (12) for light 
{C = 0) we obtain 



u 



J2 



+ 2m'U^ + + 27ti. 



(19) 



Here, a prime denotes differentiation with respect to (p, u = ^ and if is the usual angle representing 
the deflection. If we differentiate the last equation with respect to (p again we find 



fu 



+ (1 - u = 3mn^ + 7. 



(20) 



This equation with /3 = represents the deflection of light rays in Schwarzschild space-time. We also 
note that = 3-^ < is small where Rg and Rc are the Schwarzschild and compact object radii 
respectively. For example, for a compact object like the Sun this is of the order of 10""^. We therefore 
have, retaining only the first term on the right hand side of equation (20) 



3m 



Ul 



2R^ (1 - /32) 



1 + - cos 2 f 1 
3 V 



1/2 



(21) 



Equation (20) has also an exact solution if we only retain the second term on the right hand side. 
Therefore, the general solution is given by 



^ = Y3;g2+^sin(l 



3m 



2i?2 (1 - /32) 



l + ^cos2(l 



(22) 
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where R is the shortest distance between the object and the undeflected path of the hght rays in flat 
space 

R = rsmip. (23) 
At large distances, r — ^ oo, (p — ^ (^oo ^ Ij £^nd from equation (22) we obtain 

- ^ 

(1 - /^2)3/2 

As can be seen, the total deflection is twice that of the absolute value of (24) 

S = ^72^ + 27 (25) 

(1-/32)3/2 R '(l-/?2)3/2 

The first term in (25) with /3 = is the deflection angle for the Schwarzschild metric [23, 24, 33] and 
the second term is the result of the modification of the model. As can be seen, the proper potential 
increases (decreases) for the bending angle if we choose 7 > (7 < 0) respectively. 



2m 



+ 7 



(24) 



6 Radar retardation 



We now calculate the time it takes for a photon to travel from one point to another in the gravitational 
field of a central object. The equation which we shall use here is (12) with C = 0. At the point of 
the closest approach r = tq, we have dr/dt = 0. Equation (12) then yields 



(26) 



where A(r) is given by equation (7). Using equations (10), (12) and (26) we get the time taken by 
the photons traversing the distance from ro to r 



dr 



ro 



A{r) 



1 



\ r J A(ro) 



1/2- 



After substituting A{r) in the above integral, we expand the integrand to first order in 
a2r2 which are much less than unity and find 



^ / 

Jro 



This integral can be evaluated to give 
t ~ l + /32 + a2 

+ 2m In 




+ 7 



1/2 



r2 + ro^ + rro 



dr. 



(27) 
jr and 

(28) 



+ m 



r + ro 

1/2 



Vr2 — 



r + ro 



(29) 



For a = /3 = m = 0, we get the time traveled by photons in a straight line between r and ro in a 
Minkowski space-time. For a = /3 = we get the Schwarzschild metric result [23, 24, 33] and for 
7 > we find an increase in the time delay for photons compared to Minkowski and Schwarzschild 
space-times. However, in the case 7 < 0, an increase in the time delay will happen only on scales for 
which the term containing is dominant over the term including 7. 
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7 Limitations on 7 from experiment 



Let us discuss the possible range and values of the constant 7. To do this, we suppose that the value 
of c? is much less than 7 and hence the term including o? will appear only at very large scales such as 
galactic and cosmological scales. To determine the range of 7 we consider some of the tests performed 
in our Solar System. 



7.1 Deflection of light rays in the Solar System 

In Solar System tests, the Sun is the central object and wc neglect the effects of other planets. The 
second term in equation (25), = ^^2^^/^ ~ 27i?, is a simple modification to the standard result, 

namely ^qb. = 4 ^^^ with G, Mq being the gravitational constant and the Sun's mass respectively 
and c is the speed of light. The deflection of a ray that comes from infinity and grazes the Sun's 
limb is 6gr ~ 1-75 arcsec. To date, the best measurements on the deflection of light from the Sun 
have been obtained using very-long-baseline interferometry data, measuring the deflection of photons 
emanated from distant compact radio sources. Since no deviation from general relativity has been 
reported [25], we find that \S^\ < 2.11 x 10~^° and consequently I7I < 1.52 x 10~^^cm~^. Such 
experiments constrain the value of 7, but give no clue as to the sign. 



7.2 Time delay of signals in the Solar System 

According to equation (29), photons are delayed by the curvature of space-time characterized by the 
line element (6). To determine the range of 7, in this section, we compare the experimental data of 
time dilation of signals with predictions afforded by theory. In 1979 [26] the result of Viking Relativity 
Experiment confirmed the "Shapiro" time delay in the Solar System to an accuracy of 0.1%, but the 
most recent result is the frequency shift of radio signals to and from the Cassini spacecraft as they 
passed near the Sun [27] . In general relativity, the increase in At produced by the gravitational field 
of the Sun over the time taken for photons to travel the round trip between the ground antenna and 
the spacecraft at distances Vg and respectively from the Sun, is [28] 
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(30) 



where G is the gravitational constant, Mq is gravitational mass of the Sun and b is the impact 
parameter. It is convenient to use the relative change in the frequency which is caused by the 
gravitational time delay [29], because the Doppler shift due to the receiver's motion has no effect 
owing to the cancelation at both the receipt and emission of the radio signals [29]. This frequency 
shift is defined as y = — Indeed, the frequency shift was used by the Cassini experiment. For 
a case of 5 <C re,rs, which is valid for the Cassini experiment, the general relativistic contribution is 
expressed as [28] 



Mq db 



(31) 



We take the extra term caused by 7 in time delay in equation (29) and find the extra frequency shift 
as 



^7 = 7 



ro + 2re 



+ 



ro + 2r, 



(ro -I- re)2 (ro + r^) 



db 



(32) 



where we have assumed ^ ^ ^ ~ ^ near the solar conjunction (6 <C re,rs). For a spacecraft 
much farther away from the Sun than the Earth, ^ is not very different from the velocity of Earth 
Ve = ?>Qkms~^ . Now we can put constraint on 7 from the Doppler tracking of the Cassini spacecraft 
while it was on its way to Saturn reported in [27]. From [27] we find that |y-y| < 10~^^ and therefore. 
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at rs = 8A3AU, 

bmin — I-S-Rq we obtain I7I < 10 "^^cm ^. We gain seven order of magnitude more 
accuracy than the constraint from the bending of light experiments. However, it is still not possible 
to determine the sign of 7. 

7.3 Precession of planetary orbits 

Now, the measured perihelion shift of Mercury is known accurately; after the perturbing effects of 
other planets have been accounted for, the excess shift is known to be about 0.1 percent from radar 
observations of Mercury between 1966 and 1990 [30, 31]. The prediction of general relativity for 
perihelion shift of Mercury is AtpoR = 42.98 arcsecj century [23, 24] while the observed precession 
of the perihelion of Mercury is Acpobs = 43.13 it 0.14 arcsec/ century [32]. Therefore the difference 
Sip = Aifobs — AipGR = 0.15 arcsec/ century can be used to constrain 7. Prom equation (18) we 
obtain the contribution of 7 term to the perihelion shift 

where we have assumed that 3m <^ and ^ = a(l — e^), with a being the semi-major axis and the 
eccentricity of planet respectively. By assuming that the difference 6ip is due to the contribution of 
the 7 term in metric (6) (note that is much smaller than 7), the observational result imposes the 
following constraint on 7 

Use of the observational data for Mercury, equation (34), then gives I7I < 1.33 x 10"^*^ cm~^. This 
result is two order of magnitude smaller than that given by time delay of signals. 

7.4 Mean motion 

Due to the extra terms 27r and a^r^ in the metric, the radial motion of a test body around a central 
object will be affected by an additional acceleration. Let us consider a circular orbit 



ruj = — 2 7c — a c r, (35) 



GM 
r 

where lo is the angular frequency of the orbit. We can rewrite the last equation as 



Here Me// is the effective mass of the central object which can affect the motion of the test body. 
Obviously, a positive 7 would decrease the mass of the central object, leading to an excess in the 
orbital semi-major axis of the test body. In other words, the mean motion of the test body n = ^J'^- 
is changed by [36] 



— = -;^(aV + 7r2), (37) 
n 2rg ^ ' 

where r, = ^ is the gravitational radius of the central object. 

We can evaluate the statistical error on the semi-major axis of each test body 8a = — |a^, and 
interpret it as the uncertainty in the determination of 7. Bounds on 7 are shown in Table 1 where 
ba is the statistical error in the orbital semi-major axis of each planet [37]. We see that the value 
of 7 is between 3 to 5 order of magnitude smaller than that given by the time delay of signals and 
about three order of magnitude smaller than that of the perihelion shift, but since we are considering 
a statistical error we have in effect constrained the absolute value of 7. 
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oa(cm) 


{Vlim{cm ^) 


M ercury 


1U.5 


o oo v> 1A — 3'2 


Venus 




1 1 cc 1 n— 32 
i.ib X iU 


IT' J. I. 

hartn 


14. D 


1.93 X 10 


Mars 


00.7 


2.46 X 10 


Jupiter 


Doyu.u 


D.UU X iU 


Saturn 


4222.0 


6.35 X 10-32 


U ranus 


38484.0 


7.15 X lO-^"-^ 


Neptune 


478532.0 


2.32 X 10"^^ 


Pluto 


3463309.0 


7.44 X 10"^^ 



Table 1: Limits on 7 due to anomalous mean motion of planets in the Solar System. (7)am is the 
extreme limit on the parameter 7. 



8 Conclusions 

In this paper we have studied the classical tests of general relativity in a brane-world scenario in 
which the localization of matter on the brane is through the action of a confining potential. We have 
studied the effects of and 7 parameters on radar retardation, bending of light rays and the change 
in the mean motion of planets and test bodies. These are similar to the familiar results obtained 
in a Schwarzschild space-time but with certain corrections. To determine the value of constant 7, 
we assumed that the parameter is much less than 7 and hence the terms including could be 
neglected. Using available data from the Solar System, we obtained different constraints on 7. The 
best constraints come from anomalies in the mean motion of the Solar System planets resulting in 
absolute values for 7 of 1.93 x 10-^^ cm"^ and 2.46 x 10-^^ cm"^ for the Earth and Mars respectively. 
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